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Nonequilibrium molecular dynamics simulations are used 
to study the shear thinning behavior of immiscible symmetric 
polymer blends. The phase separated polymers are subjected 
to a simple shear flow imposed by moving a wall parallel to 
the fluid-fluid interface. The viscosity begins to shear thin at 
much lower rates in the bulk than at the interface. The entire 
shear rate dependence of the interfacial viscosity is consistent 
with a shorter effective chain length s* that also describes the 
width of the interface. This s* is independent of chain length 
N and is a function only of the degree of immiscibility of the 
two polymers. Changes in polymer conformation are stud- 
ied as a function of position and shear rate. Shear thinning 
correlates more closely with a decrease in the component of 
the radius of gyration along the velocity gradient than with 
elongation along the flow. At the interface, this contraction 
of chains is independent of N and consistent with the bulk 
behavior for chains of length s* . The distribution of confor- 
mational changes along chains is also studied. Central regions 
begin to stretch at a shear rate that decreases with increas- 
ing N, while shear induced changes at the ends of chains are 
independent of N. 

I. INTRODUCTION 

The viscosity of a polymer melt decreases when the 
melt is subjected to a sufficiently large shear rate, a phe- 
nomenon known as shear thinning. This non-Newtonian 
behavior is important in polymer processing and applica- 
tions and has been extensively studied both experimen- 
tally and theoretically The basic origin of shear 
thinning is that polymers elongate and align with the 
flow when they are sheared more rapidly than they can 
relax. As the shear rate increases, the degree of align- 
ment rises and the polymers present a decreasing hydro- 
dynamic resistance to flow. Recent numerical simulations 
have investigated a variety of issues in this non- 
Newtonian regime, including the shear-rate dependence 
of the viscosity, polymer conformation, and wall slip. 

One area which has yet to be explored is shear-thinning 
near the interface between immiscible polymers. Blend- 
ing of polymers is common in industry, as this is one way 
of creating new materials. However, the blended poly- 
mers are often immiscible, leading to phase separation. 
The resulting interfaces can have a major impact on the 
processing and ultimate properties of the blend. 

The static properties of interfaces in binary poly- 
mer blends are typically described by the Flory-Huggins 



model H . The degree of immiscibility is characterized by 
a parameter x, which represents the free energy cost for 
placing a monomer of one type into a homogeneous re- 
gion of the other type. There is a characteristic length s* 
of the loops of one type of polymer that penetrate into 
the region occupied by the other polymer. By balanc- 
ing the enthalpic cost and entropy gain, one finds that 
s* ~ 1/x JrJ- The interface width scales as the radius 
of gyration of a segment of length s* . 

De Gennes and coworkers argued that the dynamic 
properties of the interface should also be controlled by 
s* rather than the chain length N |]l0|. In particular, 
they suggested that for polymers in the Rouse limit the 
interfacial viscosity r// should be given by the bulk vis- 
cosity of chains of length s* . This implies that rjj is only 
a function of \ and has no dependence on N (as long it is 
larger than s*). Goveas and Fredrickson reached similar 
conclusions from more detailed calculations Q . 

In a recent paper we presented a simulation study of 
the interfacial viscosity in the Newtonian regime of a 
symmetric binary blend of Rouse polymers . As pre- 
dicted (l^Jl^], we found that the interfacial viscosity is 
determined by the degree of immiscibility and converges 
to a chain length independent value at large N. More- 
over, both t]i and the interface width decrease with in- 
creasing immiscibility and their values are consistent with 
a common value of s*. 

In this paper we investigate the shear thinning behav- 
ior of the interfacial viscosity. We find that the same s* 
that sets the width of the interface and the Newtonian 
viscosity determines the shear rate where shear thinning 
begins. Indeed the entire variation of rjj with shear rate 
follows the bulk shear thinning curve for chains of length 
s* and is independent of N. The shear rate dependent 
viscosity converges to the bulk behavior for chains of 
length A" at distances of order R g from the interface. 
Since N is larger than s*, the bulk viscosity 7]b is higher 
than rjj in the Newtonian regime and begins to shear thin 
at a lower shear rate. Indeed, the interface does not be- 
gin to shear thin until rjs has decreased to a value that 
is comparable to the Newtonian limit of r/j. At higher 
shear rates the viscosity is nearly independent of distance 
from the interface. 

We also examine changes in the conformational prop- 
erties of chains near the interface and in the bulk as a 
function of shear rate. In equilibrium, the bulk values of 
the mean-square components of the end-to-end vector are 
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equal, R 2 . = R 2 = R 2 . The two components in the plane 
of the interface (x and y) remain nearly unchanged at the 
interface, but the normal component decreases because 
of the constraint imposed by immiscibility. As the shear 
rate is increased, polymers elongate along the direction 
of flow (x), and contract in the orthogonal directions. 
We find substantial elongation before there is any mea- 
surable change in bulk viscosity. The onset of bulk shear 
thinning appears to be more closely correlated with the 
onset of contraction in the z direction. The correlation 
between viscosity and contraction is even more dramatic 
at the interface. Like plots of R 2 Z as a function of 
shear rate are independent of N and follow the bulk be- 
havior for chains of length s* . In contrast, changes in 
elongation depend on ./V and begin at lower shear rates 
than changes in 

The remainder of the paper is organized as follows. 
In the next section we review the interaction potentials 
and simulation techniques used in this work. Results and 
analyses are described in Sec. njl and a summary and 
discussion are presented in Sec. pV|. 



II. MODEL 

The polymer model and simulation techniques are sim- 
ilar to those used in previous work [p^I4| . The polymer 
potential is based on the bead-spring model developed 
by Kremer and Grest jig) . Linear polymers containing 
N beads each are created by linking nearest neighbors on 
a chain with the potential 
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where r%j is the distance between beads i and j, Rq = 
1.5a, k = 30e/cr 2 , and a and e set the length and energy 
scales, respectively. All particles in the system interact 
through a truncated Lennard- Jones potential 
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where the interaction energy t a p and cutoff r c a g depend 
on the types a and (3 of beads i and j, respectively. The 
cutoff is set at r c a a = 2 1 / 6 (7 to produce a purely repul- 
sive interaction between beads. For interactions between 
polymers of the same type e aa = e, with a = A or B. 
The two types of polymers are made immiscible |l6| by 
increasing the repulsive energy between unlike beads to 
(ab = £ba = e(l + e*), where e* = 1.2 or 3.2. Both 
values of e* are sufficient to produce phase separation for 
the range of chain lengths N considered here, N = 16, 
32, and 64. 

The bounds of the simulation cell are periodic in the 
x and y directions, with periods L x = 38.5cr and L y = 



33.4cr, respectively. In the z direction the cell is bounded 
by top and bottom walls. Each contains Nyy = 3200 
atoms tied to the sites of a (111) plane of an fee lat- 
tice by harmonic springs of stiffness n = 1320e/cr 2 . The 
walls are separated by L z = 47er. They confine 49152 
polymer beads, yielding a bead density of p = 0.8er~ 3 in 
regions far away from interfaces. Density oscillations are 
induced within a few a of the walls |17|| , However, the 
solid/polymer interfaces are only introduced to produce 
shear and are not of direct interest here. We thus restrict 
our discussion to regions more than five a from the walls, 
where wall-induced ordering is negligible. 

A symmetric blend is created in the following way. 
Polymers whose centers of mass are located in the bot- 
tom half of the simulation box are labeled type A and the 
remaining polymers are labeled B. The dividing plane is 
adjusted to ensure that there are equal numbers of poly- 
mers of each type. We also studied single phase systems 
with N = 10 and 14 at the same bead density. 

Shear flow is induced by moving the top wall paral- 
lel to the interface at a constant speed vw in the x di- 
rection. Due to the sharp wall/polymer interface and 
large polymer viscosity there can be a substantial dif- 
ference between the velocities of the wall atoms and ad- 
jacent polymer beads |T^,[l^,[l^-[l9| . We increased the 
wall/polymer interactions to limit this interfacial "slip". 
Atoms and beads interact through a Lennard- Jones po- 
tential with aw A — &wb = &■ For most simulations 
(-Wa = vL7 e and the cutoff is extended to r Wa = \/1-5ct 
to include part of the attractive region in the Lennard- 
Jones potential. At the highest shear rates studied, 
the wall/polymer interactions are increased by choosing 
tifA = ty/B — \/L7e, and r c = \/2I>a. 

The equations of motion are integrated using a fifth- 
order predictor-corrector method pOfl , with a time step 
St = 0.0075r, where r = a^Jm/e is the basic unit of 
time, and m is the mass of a monomer. A constant tem- 
perature of ksT = Lie is maintained with a Langevin 
thermostat |ljj]. To ensure that this thermostat does 
not bias the shear profile, the Gaussian white noise and 
damping terms are only added to the equations of motion 
for the velocity components normal to the mean flow (y 
and z) (l^j2^] . There is a characteristic time required for 
energy to flow from velocity fluctuations along the flow 
direction to the thermostatted components. When the 
inverse shear rate becomes comparable to this time, the 
kinetic energy along the flow direction is no longer ther- 
mostatted effectively. This limits the maximum shear 
rate in our simulations to about 0.2r _1 . However the 
correct ensemble for higher shear rates is a matter of con- 
tinuing debate, and such strongly noncquilibrium states 
are not accessible to experiments |^1 ,^2| . 

The local shear rate, 7, of the fluid is easily calculated 
by taking the local rate of change in the x component of 
velocity, function of z, i.e. 7 = dv x /dz. This is 

done by taking slices parallel to the x — y plane, of width 
0.095cr, and averaging the velocity of the monomers 
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within these slices. The viscosity within a slice is found 
from 
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where the shear stress P xz is constant throughout the 
system in steady state. Values of shear rate and viscosity 
presented below are averaged over ten slices. 

The wall velocity is varied from v w — 0.1a /t to 8. Oct /t. 
As shown in Table |, this leads to a variation by almost 
two orders of magnitude in the bulk shear rate eval- 
uated far from any interface. Note that this shear rate 
is the same in both fluids since we consider a symmetric 
melt. Statistical fluctuations drop with the total distance 
the wall moves and increase with increasing chain lengths 
(e.g. Fig. ||). The simulations at lower shear rate were 
sheared for more than a million time steps after equi- 
libration. Half this interval was used for higher shear 
rates. 



III. RESULTS 



A. Interface Width 

In Fig.|| we show the densities of the two types of beads 
in the region near the interface for e* = 3.2 and N = 16 
and 64. For this strongly immiscible case, the densities 
change from their bulk values to zero over a few ct. More- 
over, these density profiles are nearly unaffected when N 
is changed by a factor of 4 and the bulk shear rate is 
changed by a factor of 100. 

Previous studies of equilibrium interfaces have related 
the interface width to the radius of gyration o f poly mer 
segments that enter the interfacial region 
The length of these segments, s* , is determined by the 
degree of immiscibility and becomes independent of N 
in the large N limit. Earlier studies of the equilibrium 
properties of the model considered here [jl2j are consistent 
with these predictions. For example, for e* = 3.2 the 
density profiles are independent of N for N > 16. As 
illustrated in Fig. [j], this behavior extends to the highest 
shear rates studied here. The same insensitivity to shear 
rate is found for the less immiscible case of e* = 1.2 where 
the interface is slightly wider. 



B. Viscosity 

The bulk viscosity of polymers whose length is shorter 
than the entanglement length can be described with 
Rouse theory (g5|. The limiting Newtonian viscosity at 
low shear rates scales linearly with chain length: t\b = 
C,b 2 pN, where £ is the monomeric friction coefficient and b 
is the statistical segment length. DeGennes and cowork- 
ers jl0| argued that the viscosity in the interfacial region, 
7]i, should be determined by the effective length s* rather 



than TV because s* is the length of segments that must 
relax during shear flow. Their conclusion is supported by 
the analysis of Goveas and Fredrickson pd| and by our 
earlier simulations in the Newtonian limit]l2| . 

The dependence of the bulk and interfacial viscosities 
on chain length in the Newtonian limit is illustrated in 
Figure ||. Here 77 is plotted against height z for e* = 
3.2 over the range of chain lengths where the interface 
width, and thus s*, is constant. Away from the interface 
77 approaches the bulk viscosity, which rises linearly with 
N for N up to 64. The slope of this rise determines 
C = 0.38±0.03to/t, since b — 1.28ct is known from static 
properties [JL2 15 and p = 0.8ct~ 3 . 



In the interfacial region the viscosity is smaller than 
the bulk value and is independent of chain length. This 
is consistent with the picture that the constant value of 
s* determines rjj. Equating s* to the chain length that 
would give a bulk viscosity equal to r\i yields s* ~ 10. 
A similar analysis for e* = 1.2 yields s* w 14, and in 
both cases the radius of gyration corresponding to s* is 
1.6 times the interface width [[l^ll . Thus the Newtonian 
response is consistent with the picture advanced by de- 
Gennes and coworkers Jl0| . 

The shear thinning behavior of the bulk viscosity with 
increasing shear rate is shown in Fig. ||. As is well known, 
the onset of shear thinning occurs when the system is 
sheared more rapidly than it can relax. Since the time 
to relax increases with chain length, the onset of shear 
thinning moves to lower shear rates with increasing N 
0. As in previous work (7J, the decrease in viscosity 
at high shear rates can be fit to a power law 77 ~ 7~ a , 
with a near 1/2. The dashed line in Fig. || shows a 
fit to N = 64 data with a — 0.47. The low shear rate 
viscosity increases with chain length, and the high shear 
rate behavior is nearly independent of N . 

The variation of the interfacial viscosity with the in- 
terfacial shear rate is shown in Fig. |i| for different N 
and e*. In all cases, r\i begins to shear thin at much 
higher shear rates than rjg. Data is only presented for 
values of N where the Newtonian response and interface 
width are independent of chain length, so that s* should 
also be constant. One sees from the figure that the shear 
thinning behavior is also independent of N in this limit. 
Moreover, the entire shear thinning curve for all N at 
each value of e* is consistent with the bulk shear thinning 
behavior of chains with length equal to the values of s* in- 
ferred from the corresponding Newtonian response. This 
is strong evidence that a single segment length deter- 
mines the non-Newtonian shear thinning behavior as well 
as the static interface profile and Newtonian response. 



C. Polymer Conformations 

We now describe the changes in chain conformation 
that are associated with the shear thinning of the bulk 
and interface. Fig. ^| shows the mean square of the three 
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components of the end-to-end vector J? as a function of 
the height of the polymer's center of mass. The two pan- 
els show results for two low shear rates with N = 16 
and e* = 3.2. The lowest shear rate (panel (a)) is well 
into the Newtonian regime, and the polymers have time 
to relax to equilibrium conformations. Far from the in- 
terface the polymers follow isotropic random walks with 
R 2 = Ry = R 2 w 8fi 2 . Near the interface R 2 is strongly 
suppressed because immiscibility eliminates paths that 
take polymers too far into the other phase. The values of 
R 2 and Ry are nearly unaffected by this constraint, but 
show a small rise that is comparable to our statistical 
noise. The value of R 2 does not attain its bulk level until 
the center of mass is far enough from the interface that 
chains are unlikely to reach it. This distance is compa- 
rable to the bulk value of \JR%, which is larger than the 
range over which the density changes (Fig. yj). 

Panel (b) of Fig. [| shows the mean square components 
of R at the third lowest shear rate in Fig. [| Surprisingly, 
although the viscosity is nearly indistinguishable from 
the low shear rate limit, there is a substantial change in 
conformation. The component along the flow direction, 
R 2 , increases by about 25% throughout the system. The 
bulk value of R 2 Z drops by about 10%, while the interfacial 
value is nearly unchanged. This is consistent with the 
observation that the interface shear thins at higher shear 
rates than the bulk. 

Conformation changes at much higher shear rates are 
illustrated in Fig. [| Here the mean square x and z com- 
ponents of the end-to-end vector are plotted as a function 
of center of mass height for N — 32, e* — 3.2, and four 
values of vw- As expected, chains in the bulk regions 
are stretched and aligned along the flow direction. At 
the lowest shear rate, the bulk value of R 2 is about 20% 
above its equilibrium value. At the highest shear rate, 
the length of the polymer along the flow direction has 
grown to about 40% of the fully extended length. The 
components orthogonal to the flow decrease. The drop 
in R 2 is much larger than that in Ry, which drops from 
about 17 to 10<7 2 over the studied range of 7. Overall 
there is a net increase in R 2 from about 50 to 180<r 2 , be- 
cause the increase in R 2 is much larger than the decrease 
in the other two components. 

Chains in the interfacial region also stretch along the 
flow direction and contract in the orthogonal directions. 
However, as with the viscosity, the shear rate dependence 
is shifted to higher 7 at the interface. The value of R 2 
is initially slightly higher at the interface and increases 
more slowly with 7, leading to a pronounced dip at the 
highest shear rates. In contrast, the dip in R 2 Z at the 
interface decreases with increasing shear rate because the 
bulk value drops more rapidly with 7. 

Fig. |^ provides a more detailed picture of the variation 
in the conformation of chains near the interface with in- 
creasing shear rate. Interfacial values of R 2 Z are plotted 
against shear rate for N = 16, 32, and 64 at e* = 3.2. 
Bulk values for chains with length equal to the inferred 



value of s* = 10 are shown for comparison. All curves are 
flat at low shear rates, and begin to decrease at the point 
where the interfacial viscosity begins to fall in Fig. || As 
with the data for 777 , the interfacial conformations for all 
chain lengths collapse on to a universal curve that coin- 
cides with the bulk curve for chains of length s* within 
our errorbars. This provides further evidence that a sin- 
gle time scale corresponding to s* controls the behavior 
of the interface. Note that the uncertainties are larger 
for interfacial values due to the relatively low number of 
polymers whose center of mass lies at the interface. As 
discussed previously |l^,^], there is a preponderance of 
chain ends at the interface. This reduces the number of 
chain centers at the interface by a factor of three for the 
parameters considered here. 

The idea that relaxation times should scale with the 
length of a subset of the chain rather than its entire 
length is related to other recent work . These papers 
have examined changes in the conformation of individual 
chains pulled through gels or other polymers. Segments 
at the free end of the chain are assumed to relax with 
a time scale that is characteristic of bulk polymers with 
the same length as the segment. Since the relaxation time 
grows with the length of the segment, the polymer is less 
relaxed and more stretched as one moves away from the 
free end toward the pulling end. 

A similar phenomenon should occur within individual 
chains in a sheared polymer melt. Segments at the end 
of the chain will be more relaxed than the center leading 
to a dumbbell configuration that narrows in the center. 
In addition, the ends should begin to stretch at higher 
shear rates than the center because their characteristic 
relaxation time is shorter. Fig. ^ shows the end-to-end 
distance of segments of length 4 as a function of shear 
rate. Results for the two end-segments are averaged and 
compared to a segment in the middle of the chain. Note 
that the ends show the same behavior for all ./V indi- 
cating that the characteristic relaxation time for their 
conformations is only a function of segment length. In 
contrast, the middle segment becomes more stretched as 
N increases and the characteristic rate at which it can no 
longer relax decreases with increasing N. These effects 
are reminiscent of a nematic coupling effect seen recently 
in polymers, where the orientation of middle segments 
also relaxed more slowly with increasing N, while the 
rapid relaxation of ends was independent of N p8] |. 

IV. SUMMARY AND DISCUSSION 

In this paper we have examined shear thinning and 
conformational changes in the bulk and interfacial re- 
gions of a phase-separated binary blend of Rouse chains. 
The bulk behavior is consistent with previous studies. 
The most surprising result is that pronounced elonga- 
tion of the chains is observed before there is a notice- 
able change in viscosity. The onset of significant changes 
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in t]b seems more closely connected to the contraction 
along the gradient direction, which becomes appreciable 
at higher shear rates than the elongation along the flow 
direction. 

The changes in intcrfacial viscosity are consistent with 
the prediction |l(],[llj that the interfacial viscosity is de- 
termined entirely by the characteristic length s* of loops 
that cross the interface between immiscible polymers. 
This length is determined by the degree of immiscibil- 
ity and independent of N for the polymers studied here. 
The interfacial viscosity shows these same trends. For 
each degree of immiscibility, plots of interfacial viscos- 
ity against interfacial shear rate for different N collapse 
onto a single curve. The entire curve coincides with the 
bulk shear thinning behavior of polymers of length 10 for 
e* = 3.2 and 14 for e* = 1.2. Thus a single effective chain 
length describes both the Newtonian value of r\i and its 
shear thinning behavior. Our previous work showed that 
the same effective chain length described the width of 
equilibrium interfaces. Jl2[ 

The suppression of the viscosity near the interface (Fig. 
|J) leads to an effective slip boundary condition [jni] p^[ . 
The amount of slip is most pronounced in the Newtonian 
limit where the difference between r\i and t\b is largest. 
As the bulk viscosity begins to shear thin, the amount 
of slip decreases. Once the shear rate is high enough to 
produce significant shear thinning at the interface, the 
bulk and interfacial viscosities converge (Fig. [|). In this 
limit the amount of slip is negligible. 

Changes in the conformation of polymers near the in- 
terface were also studied. In equilibrium, the interface 
suppresses the component of the end-to-end vector per- 
pendicular to the interface, but has little effect on the 
in-plane components. The changes in interfacial viscos- 
ity are most closely correlated with changes in the per- 
pendicular component. Like 77/, plots of R\ against shear 
rate for different chain lengths collapse onto a common 
curve that is consistent with that for bulk chains of length 
s* (Fig. 0). While the elongation increases more slowly 
near the interface than in the bulk, it begins at much 
lower shear rates than the change in interfacial viscos- 
ity (i.e. Fig. ||) an d depends on chain length. Thus 
our results for both the bulk and interface indicate that 
polymer contraction perpendicular to the interface is the 
most important structural change associated with shear 
thinning. 

We also explored changes in the conformation of seg- 
ments of 4 beads within longer chains as a function of 
shear rate. Segments in the center of the chains began to 
elongate at the same shear rate where the bulk viscosity 
for the given chain length showed shear thinning. How- 
ever, the conformation of the ends remained unchanged 
until much higher shear rates. Moreover, the conforma- 
tion of the ends was nearly independent of the total chain 
length and their shear thinning behavior was comparable 
to that for bulk chains with length 4. Recent studies of 
much longer chains have used the concept of a position 
dependent relaxation time to determine the conformation 



of a single chain in a solvent or gel |27| . Our study shows 
that the same concept applies to polymer melts. 
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TABLE I. The shear rate in bulk regions of each polymer 
7s at the indicated values of chain length and wall velocity 
vw ■ Uncertainties are less than 5%. The proportionately 
larger shear rates at vw = 8a /t are due to an increase in the 
wall-coupling parameters, as discussed in Sec. 



Vw 






7s 








TV = 10 


TV = 14 


TV = 16 


TV = 32 


TV = 64 


0.1 


.0019 


.00175 


.00159 


0.0013 


0.00094 


0.3 


.0055 


.00512 


.00479 


0.0039 


0.0038 


0.5 


.0091 


.00872 


0.0081 


0.0074 


0.0073 


1.0 


.018 
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FIG. 1. Densities of type A (decreasing curves) and B (in- 
creasing curves) beads showing that the interface width is 
insensitive to chain length and shear rate. The immiscibility 
parameter is e* = 3.2 and circles are for TV — 64, vw = 0.1; 
diamonds are for TV = 64, vw = 5.0; and triangles are for 
TV = 16, vw = 5.0. Corresponding values of 7b are given in 
Table @. 




FIG. 2. Newtonian viscosity as a function of height z for 
e* = 3.2. The bulk viscosity scales linearly with chain length 
for TV = 16 through 64, while the interfacial viscosity r\i is 
independent of chain length over this range. 
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FIG. 3. Bulk viscosity r\B as a function of shear rate 73, 
for polymers of length N = 16, 32 and 64. As N increases, 
shear thinning is observed at lower values of 7s. The dashed 
line shows a power law fit t\b ~ j~ ' 47 . 
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FIG. 4. The interfacial viscosity r/j as a function of inter- 
facial shear rate 71 for (a) e* — 1.2 and iV = 32 and 64, and 
(b) e* = 3.2 and N = 16, 32 and 64. The bulk viscosities 
of polymers with length equal to the values of s* determined 
from the Newtonian response are shown for comparison. The 
dashed lines show the bulk viscosity of chains of length 64. 
Interfacial shear rates are evaluated over an interval of width 
one a in the center of the interface. 




FIG. 5. The mean square components of the end-to-end 
distance along a = x, y, and 2 as a function of the cen- 
ter of mass height z for N = 16 polymers with e* = 3.2 at 
(a) vw = 0.1<t/t and (b) vw = 0.5<t/t. Although the vis- 
cosities at both shear rates are nearly indistinguishable from 
the limiting Newtonian viscosity, there is a clear stretching of 
polymers along the flow direction in (b). 
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FIG. 6. The mean square x and z components of the 
end-to-end vector as a function of center of mass height for 
N — 32 polymers at wall velocities vw = 0.1a /t (diamonds), 
0.5a /t (circles), 2a /r (squares) and 5a /t (triangles). Corre- 
sponding bulk shear rates are given in Table |l| 
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FIG. 8. The mean square end-to-end distance of 4 bead 
segments at the middle (filled symbols) and ends (open sym- 
bols) of polymers with the indicated TV as a function of shear 
rate. The ends show the same behavior for all TV. The mid- 
dles begin to stretch at a lower shear rate that decreases with 
increasing TV. 
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FIG. 7. The mean square z component of the end-to-end 
vector, R^, for polymers centered at the interface with 
TV = 16, 32 and 64. Bulk values for chains with TV = 10 
are shown for comparison. All curves begin to decrease at the 
same shear rate, and this shear rate coincides with the drop 
in interfacial viscosity. 
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